Modification of anatase TiO2(001) surface electronic structure by Au impurity 
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We have used density functional theory calculations based on the projector augmented wave 
method to investigate the electronic structure of Au-incorporated anatase TiO2(001) surface. Due 
to the coordination with several level oxygens, Au atoms can be encapsulated inside Ti02 slab. 
Au is adsorbed over the surface Ti-O bond, so called the bridge site on anatase TiO2(001)-lxl 
surface. However, for 0.25 ML coverage, Au atoms energetically prefer to stay at 0.64 A above 
the midpoint of the two surface oxygens which is significantly closer to the surface layer. When 
implanted inside the slab for full coverage, Au forms parallel metallic wires inside Ti02 lattice 
where interlayer distances increase due to local segregation. Au brings half-filled impurity states 
into the band gap leading to metallization, in addition to other filled surface and impurity bands 
within the gap. These Au-driven Fermi-level-pinning gap states are close to, or even in some cases 
inside, the conduction band of the host slab. On the other hand, if Au is substituted for the surface 
Ti atom, Fermi level falls lower in the gap closer to the valence band top. 

PACS numbers: 73.20.Hb, 68.43.Bc 



I. INTRODUCTION 

Lately, titania (T1O2) has received an increased at- 
tention since it is considered to be promising for cost- 
effective photovoltaic applications. However, high reac- 
tivity of TiC>2 under only UV light bears a great dis- 
advantage which reduces the quantum efficiency so that 
pure Ti02 is not alone sufficient for a practical system 
applicationfii£iMi£ 

McFarland and Tang£, in a recent work, proposed 
an Au/Ti02/Ti multilayer photovoltaic device on which 
photon absorption occurs in the deposited dye molecules 
while electron-hole pairs are created upon light illumina- 
tion inside the semiconductor for the conventional solid- 
state solar cells. This ultrathin metal-semiconductor 
junction Schottky diode has driven a particular atten- 
tion to gold-titania interface^. 

Gold has been proposed to enhance the catalytic ac- 
tivity of anatase TiO^*^ Au incorporated anatase-bascd 
nanocatalysts have been synthesized for device appli- 
cationSi 10 ! 11 ' 12 Moreover, diffusion of gold into anatase 
polymorph of titania has also been reported^ 

The anatase phase of titanium dioxide, although be- 
ing less stable than the rutile polymorph, is catalytically 
more activ e) 14 ' 15 : 16 . In this sense, the surface properties 
are of major importance. There are several studies on the 
topological and electronic structure of the single crys- 
tals of anatase TiO 2 (001) surfaceiii^ 1 ^. In this pa- 
per, we investigated the reconstructive effect of various 
types of Au incorporation in the anatase substrate and 
the role of such an impurity on the electronic structure 
of TiO2(001). We considered Au in and on the surface 
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at quarter and full monolayer (ML) concentrations in- 
cluding Au(Ti) substitutional cases as well as gold dimcr 
adsorption on the surface. We have discussed their ther- 
modynamic stabilities for various experimental environ- 
ments. 



II. METHOD 

The total energy density functional theory (DFT) 
calculations have been performed within the general- 
ized gradient approximation (GGA) for the exchange- 
correlation effects via Perdew-Burkc-Ernzcrhof (PBE96) 
functional 21 -, using plane-wave basis sets and the projec- 
tor augmented waves (PAW) metho d 22 ' 23 as implemented 
in the Vienna Ab-initio Simulation Package (VASP) 2 ^. 

Anatase TiO2(001) stoichiometric surface has been 
modeled as a symmetrical slab which constitutes 6 Ti02 
layers in a supercell with a vacuum region of at least 13 A. 
A Ti02 layer consists of 3 atomic planes in which bridging 
oxygen atoms are out of the level Ti-plane. Our choice re- 
garding the number of layers represent one of the thickest 
slab models among the other theoretical studie s 17 ' 18 : 19 : 20 . 
Moreover, for Au implantation deeper than the fourth 
atomic plane we used 8-Ti02-layer slab model to avoid 
any interaction between the Au impurities within the 
same supercell which are implanted from the bottom and 
from the top surfaces. 

For the gold incorporated surface models we kept ions, 
only in the central Ti02 layer, fixed to their bulk posi- 
tions because the forces on these remain to be insignif- 
icantly small in all of the calculations. We employed a 
conjugate-gradients algorithm to compute the electronic 
ground state based on the minimization of the total en- 
ergy subject to a convergence tolerance of 0.1 meV. Ge- 
ometry optimization has been performed, subsequently, 



2 



by the reduction of the quantum force on each of the 
unconstrained ions to less than 10 me v/A. 

Single particle wavefunctions have been expanded in 
terms of plane waves up to a cut-off energy value of 400 
cV. The Brillouin zone (BZ) integrations have been car- 
ried out over 32 and 5 special fc-points sampled in the 
irreducible wedge for lxl and 2x2 surfaces, respectively. 
We did not impose any symmetry in our surface calcula- 
tions. 

In thermodynamic equilibrium, the most stable sur- 
face composition, that is entirely surrounded by a 
thermal bath at temperature T under a given pres- 
sure p, minimizes the surface Gibbs free energy j(T,p) 
(GFE),2£ 1 2S 1 212§. Therefore, relative stability of the 
phases can be expressed in terms of the difference 
A"f(T,p) in the GFE of impure surface structures and 
relaxed clean surface, as, 

A7(T,p) = — |G Au/Ti o 2 (T,p, An Ti , An Au ) - G T io 2 (T,p) 
-An T i"Ti(T,p) - AnAuAiAu(T,p)| , 



where G 



Au/TiO. 



and Gtioo are the GFEs of the Au incor- 



porated and reference clean Ti02 surface slabs, respec- 
tively. A is the corresponding unit cell area, also func- 
tioning as normalization for different stoichiometrics. As 
the name suggests, An is the difference in the number of 
atoms from that of the reference surface while /i denotes 
the chemical potential for the referring species. A nega- 
tive A7 value indicates a more stable structure relative 
to the clean Ti02 surface. Positive values, on the other 
hand, give the relative formation energy that is required 
to assemble the corresponding composition. 

GFE, in general, is defined by, G = U+pV-TS. One 
can omit pV term in comparison to the surface energy for 
the pressures, p, under consideration as well as the TS 
term since the contribution from the entropy, S, is negli- 
gible. Therefore, it can be approximately expressed only 
by the internal energy U. If one also neglects the ionic 
vibrational contributions, U will be equal to the total en- 
ergy .Etotai(-W, V) obtained by the DFT slab calculations. 

Thermodynamic equilibrium is established when the 
chemical potential of a given atomic species becomes 
equal in all phases that come into contact with each 
other. In particular, we assume that the surface struc- 
tures are in equilibrium with the anatase Ti02 bulk so 
that fj,Ti + 2/xo = A*Ti0 2 ■ Moreover, the reference ener- 
gies at the most stable elemental phases of Ti and Au set 
the upper bounds for their corresponding chemical poten- 
tials. For instance, the maximum value of uxi can be ac- 
cessed in the hep Ti bulk solid phase. Otherwise, the bulk 
phase would be unstable with respect to precipitation of 
bulk Ti. Similarly, fiAu can not be above the chemical 
potential of ccp gold bulk. Molecular oxygen defines the 
upper boundary for fiQ so that fiQ = ^Eo 2 which cor- 
responds to an O-rich experimental environment, there- 
fore, defining the minimum value of /ixi through the ther- 
modynamic equilibrium condition, fxn + 2/iq = "Ti0 2 ■ 



III. RESULTS AND DISCUSSION 

We systematically studied Au impurities on and inside 
the anatase TiO2(001) stoichiometric surface. The role 
of such an incorporation on the electronic properties of 
anatase surface has been investigated by examining ad- 
sorptional, substitutional, and interstitial impurities at 
quarter and full coverages. For these concentrations, we 
considered gold ions in the subsurface at the interstitial 
cavities up to a depth of the fifth atomic plane, and also 
substituted them for Ti ions up to the third Ti plane. 

The topological structure and the electronic properties 
of the bare anatase TiO^OOl) has been discussed in de- 
tail elsewhere^. Fig. [1] shows the optimized geometries 
and corresponding electronic structures of Au-Ti02 com- 
bined systems for lxl surface. The first two panels to 
the left show the side views of the geometric structures 
along [100] and [010] directions. In addition, relevant 
electronic structures are presented next to these geome- 
tries at each row representing a distinct impurity case. 




m x' r M 



FIG. 1: (color online) The geometric and electronic structures 
of a-1 x 1 (first row), si— 1 x 1 (second row), and s2-l x 1 (last 
row), for the Au/TiO2(001)-lxl systems. Big lightgray (yel- 
low) balls represent gold. Small dark (red) and bright (gray) 
balls denote oxygen and titanium atoms, respectively. Bond 
angles and bond lengths are given in degrees and angstroms, 
respectively. Energy bands and DOS structures are described 
in detail in the corresponding subsections. 
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Each ion is labeled with respect to the atomic plane, that 
it belongs to, for each of species separately. For instance, 
an oxygen at the second oxygen atomic layer is labeled 
as 02. Furthermore, bulk termination gives rise to un- 
dercoordinated ions by breaking the axial bonds over the 
surface. As a result, 01 and Til become twofold and 
fivefold coordinated, thus, also referred as 02c and Ti5c, 
respectively. 

A. Au/TiO 2 (001)-(lxl) 

1. Adsorptional case: a-lxl 

Au-adsorptional case, that the first row of Fig. Q] refers 
to, has been obtained by relaxing all possible initial con- 
figurations in which Au comes into contact with the ac- 
tive sites on the support surface. Truly, in this resultant 
geometry, Au relaxes to a bridge position on 01 and Til 
that reflects the maximum interaction between the Au 
adsorbate and the surface. At this minimum energy po- 
sition Ol-Au-Til angle becomes 45.4° with Au-01 and 
Au-Til bond distances being 2.39 and 2.83 A, respec- 
tively. The reconstructive effect of the adsorbate on the 
lattice remains negligibly small at the subsurface layers. 
Au adsorption at 1ML coverage adapts an ideal-like Ti02 
atomic arrangement by saturating the dangling bonds. 
Yet, the symmetry breaking in the Til-01 bonds which 
stems from the relaxation of the clean surface, is still not 
lifted as bond lengths read 2.05 and 1.86 A. The struc- 
tural parameters as well as some other physical entities 
related to electronic properties such as the surface work 
function and Fermi energy relative to bulk valence band 
maximum (VBM) are presented in Tabled] 

Au adsorbate brings a two-dimensional impurity state 
in the gap closer to the conduction band (CB). Fermi 
level passes through the saddle point at X' of this band 
corresponding to a logarithmic singularity in the site- 
projected (local) density of states (LDOS). This half- 
filled Fermi-level-pinning state, therefore, leads to metal- 
lization. The reactivity of the 01 site is higher than that 
of Til, which is also apparent from the LDOS presented 
on the rightmost panel of the corresponding electronic 
band structure in Fig [T] A group of four occupied defect 
states lie distinctly below the Fermi energy derived from 
the valence band (VB) possessing dominant 01 charac- 
ter. The contribution of Til to the total DOS around the 
Fermi energy is very weak. 

2. Substitutional case: sl-lxl 

The most noticeable effect of Au substitution for sur- 
face Ti (Til) at 1 x 1 structure, also referred as sl-1 x 1, is 
that the distance between the top and the second Ti02 
layer increases substantially. For instance, the Til-03 
axial bond length of 1.96 A at the clean surface extends 
to 2.59 A for Au-03. In addition to the increase in the 



interlayer distance, this extension also gets the Ti2-03- 
Ti2 angle to enlarge from 152.0° to 159.7° at the second 
Ti02 layer. Moreover, Au-01 equatorial bonds along 
[100] at the surface layer become almost equal being 1.95 
and 1.96 A as shown in the second row of Fig. [T] 

Unlike to that of the adsorptional case, Au-Til sub- 
stitution drives Fermi energy to fall lower in the gap 
closer to the valence band top. Because of the partially- 
occupied nature of this state which is populated by Au, 
01, and chiefly by 02, the combined system goes into 
the metallic regime. The upper lying unoccupied defect 
state, which takes up a big portion of the gap to the 
energies above the Fermi level as a result of the strong 
dispersion, shows dominant Au and 02 contribution. 



3. Substitutional case: s2-lxl 

The optimized geometry presented in Fig.[ljs2) for the 
second substitutional case (s2-lxl) has been obtained 
by relaxing an initial configuration which was formed by 
replacing Ti2 with Au in the bare 8-TiO2(001)-layer slab 
model. Au impurity sits below the surface in the second 
Ti02 layer having extended interlayer distances, similar 
to sl-lxl case, of 2.24 and 2.25 A with the top and the 
third layers, respectively. On the other hand, Til-01 
bond lengths remain to be nonequivalent with values of 
2.12 and 1.83 A making a Til-Ol-Til angle of 144.8°. 

Subsurface substitutional gold sets the Fermi level at 
0.93 cV above the valence band maximum (VBM) pinned 
by a strongly dispersed empty defect state, of ^2 eV 
wide, which descends from the CB. The second band 
below, with a width of about 0.5 eV is partially filled 
causing metallization of the system. This band is de- 
generate with a third band of similar width at T point, 
and all these three states are symmetrical with respect 
to M point along rXMXT. The unoccupied first band 
is originating mainly from the Au-04 and Au-Au inter- 
actions. The half-occupied second one shows dominant 
02 and relatively weaker 04 character. Third band has 
an LDOS peak stemming from 03. The fourth band in 
the gap which crosses the third along most of XMX' is 
due to the clean surface (see Fig. 2(a) in Rcf. [18]). 
The LDOS for the surface oxygen (01) disperses over 
VB similar to that of sl-lxl implying a charge transfer 
to subsurface cation sites. This results in a work func- 
tion of 6.06 eV. We have also considered Au substitution 
for Ti3 (not shown) . The band structure for s3-l x 1 has 
common characteristics with that of s2-lxl. The first 
and second bands are very similar, third is no longer de- 
generate with the second at T, and fourth, the surface 
state, is pushed down towards the VB, and is also sym- 
metric with respect to M. In quest of finding a trend 
we checked the band structures for Au substituted even 
deeper, namely for Ti5 (s5-lxl) and also considered a 
case for bulk substitution. Comparison of these series of 
cases, Au for Til to bulk Ti, band structures projected 
to the same surface Brillouin zone shows that the second 
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TABLE I: Calculated electronic and structural parameters for Au-TiC>2(001) anatase system: work function, position of the 
Fermi energy relative to bulk valence band maximum, Au-depth relative to surface oxygens, Au-0 and Au-Ti distances for 
each model labeled in compliance with Fig. [TJ Fig. [2] and Fig. [3] 



Model 


W(eV) 


E F (eV) 


/lAuA 


dAu — o A 




C?Au — Ti A 




a— 1 x 1 


5.45 


1.93 


2.14 


2.39(01) 




2.83(Til) 




sl-lxl 


6.37 


0.31 


-0.51 


1.95(01), 


2.03(02) 


3.51(Ti2) 




s2-l x 1 


6.06 


0.93 


-3.22 


2.24(02), 


1.98(03, 04), 2.25(05) 


3.23(Til) 


3.28(Ti3) 


a-2x2 


5.02 


2.02 


1.94 


2.20(01) 




2.68(Til) 




b-2x2 


4.29 


2.60 


0.64 


2.06(01), 


2.73(02) 


3.06(Til) 




c-2x2 


4.33 


2.59 


-2.71 


3.17(02), 


2.09(03), 2.50(04) 


2.85(Til) 


2.77(Ti2) 


d-2x2 


4.49 


2.57 


-3.40 


2.46(03), 


2.05(04), 3.13(05), 2.74(06) 


2.80(Ti2) 


2.94(Ti3) 


e-2x2 


4.47 


2.58 


-5.24 


3.11(04), 


2.06(05), 2.49(06), 3.55(07) 


2.93(Ti2) 


2.78(Ti3) 


sl-2x2 


6.80 


0.31 


-0.27 


2.01(01), 


2.03(02) 


3.37(Ti2) 




s2-2x2 


6.66 


0.58 


-3.18 


2.07(02), 


2.59(03), 2.02(04), 2.03(05) 


3.45(Til) 


3.18(Ti3) 


aa-2 x 2 


6.28 


1.04 


1.57 


2.32(Aul- 


-01), 2.07(Au2-Ol) 


2.66(Aul-Til), 2.67(Au2-Til) 



and third bands continue to repel each other, while Ef 
passes through the second defect band as before. 



^ . Interstitials at lxl 

For the interstitial case at lxl surface, strong metal- 
metal interaction distorts the local lattice structure dras- 
tically due to shortened Au-Au distance of 3.76 A. This 
causes a local segregation that leads to the formation of 
parallel metallic Au wires along [010] inside TiC>2 lattice 
where inter layer distances further increase (from ~2.02 
to 3.89 A). 



For a-2x2 case, energy gap region is characterized by 
weakly dispersing six defect states as shown in Fig. [2Ja). 
The one lying just below the CB that is half filled by 
simple electron counting drives the combined system into 
metallic state. This Au-Ol driven state is flat along JK 
and almost flat along JT and therefore quasi-one dimen- 
sional in nature as seen also in LDOS. The next defect 
state that lies 1.56 eV below the Fermi energy at T be- 
ing also quite flat, especially along JKJ', has mainly 02 
character. The remaining four states come up slightly 
above the VBM with many crossings originating mainly 
from Au-Ol and Au Til interactions. 



B. Au/TiO 2 (001)-(2x2) 

1. Adsorptional case: a-2x2 

Similar to the case for the 1ML coverage, single Au 
adsorbate at the 2x2 surface is twofold coordinated with 
01 and Til at the bridge position over the Til-01 bond 
as shown in Fig. [2^a) . Au promotes the nearest-neighbor 
Til upward elongating the Til-03 bond from a clean 
surface value of 1.96 to 2.02 A. It also interacts strongly 
with the nearest 01 which is elevated up leading to a 
significant increase in the Til-01 bondlength, from 2.16 
to 2.40 A, over which Au stays at the minimum energy 
position. As a result, Ol-Til-Ol angle gets wider by 
4° compared to the value at the clean surface. Quar- 
ter ML Au adsorption yields significant differences in the 
atomic positions from those of the 1ML coverage. For in- 
stance, Au-Til and Au-Ol bond lengths become shorter 
(2.68 and 2.20 A) than the values (2.83 and 2.39 A) for 
1ML coverage, respectively. This can be explained by 
the increased charge transfer between the impurity and 
the support surface due to much weaker Au-Au impurity 
interaction at 1/4 ML coverage. 



2. Adsorptional case: b~2x2 

In the second adsorptional model, namely b-2x2, Au 
relaxes into a very symmetrical position, 0.64 A above 
the midpoint between two 01 ions. It forms two and 
four equidistant bonds with those 01 ions and with Tils 
that read 2.06 A for the former and 3.06 A for the latter, 
respectively This isotropic coordination with two near- 
est neighbor Ols and with four surface 'Ti's get Til-03 
bonds to be aligned parallel to [001]. Hence, each of the 
Til-Ol bonds become equal to 2.10 A in length making 
and angle of Til-Ol-Til angle of 137.6°. On the other 
hand, two Ti2-02 bonds that are coplanar with Au re- 
main to be slightly skewed with congruent angles of 3.6° 
each. 

Fermi level falls into the CB by a resonant defect state 
causing metallization for the electronic structure pre- 
sented in Fig.[2{b). Two flat going occupied defect states 
lie about 0.5 and 0.9 eV above VBM. The higher one of 
these bands carries dominant 01 character, whereas the 
lower one originates from Au-Ol interaction with larger 
Au mixing, and has some 02 contribution as well. 
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FIG. 2: Geometric and electronic structures of Au/TiO2(001)-2x2 systems for the adsorptional, substitutional, and interstitial 
gold incorporations. Label for each case, is given in the DOS panel, should read as, e.g., a-2 x 2 for the first row, and so on. 
Naming convention for the atoms follow those of Fig. [1] The bond lengths and angles are given in angstroms and in degrees, 
respectively. 



3. Substitutional case: sl-2x2 

When substituted for surface titanium as shown in 
Fig. [Hsl), Au relaxes into a position, 0.38 A above the 
original Til lattice point, where it forms two Au-Ol 
bonds with slightly unequal lengths of 1.99 and 2.01 A. 
It is also twofold coordinated to two nearest-neighbor 
02s at a distance of 2.03 A away from each. The in- 
cquivalcncy in the Ti-Ol bond lengths still remains be- 
ing 1.84 and 2.06 A over the undcrcoordinated surface 
oxygen row along [100] next to the row in which Au sub- 
stituted. Moreover, Au-03 separation elongates to 2.65 
A which can be attributed to Au attaining a nominal 
charge state that weakens the charge transfer from Au to 
03. Resulting bond length of Ti2 with this 03 shortens 
to 1.85 A while the other Ti2-03 distances are all 1.96 



A. 

A stack of surface bands disperse over an energy range 
of 0.33 eV with intercrossings just above VBM. The high- 
est of those is half occupied setting the Fermi energy at 
0.31 eV relative to bulk valence band top. Therefore, 
combined system has been predicted to show metallic 
behavior. The two sharp LDOS peaks at and just be- 
low the Fermi level in Fig. [DJsl) are mainly due to 01 
which is bonded to Au. The rather lower part of those 
surface states, on the other hand, originate from under- 
coordinated Ols which are away from Au site. An empty 
defect state falls into the gap region at around 2 eV at 
T below CB having a minimum at K. This band has two 
saddle points at J and J' which are not at the same en- 
ergy because of the unequal Au-Ol bondlengths along 
[100]. Corresponding LDOS peaks come out as a result 
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of excess anti-bonding charge localized around gold. 



4- Substitutional case: s2-2x2 

Au substitution for Ti2 causes local disturbances on 
the nearby 03s as a result of the excess charge embed- 
ded by and localized around the impurity site. Equatorial 
Au-03 bonds on (100) plane extend to 2.59 A leading to 
a significant dislocation of 03s from their lattice posi- 
tions. This subsurface reconstruction induces stress on 
Til-03 interaction that flips up Til-02 bonds. Then, 
the displaced 02 levels itself with 01 atomic plane as 
shown in Fig. [2js2). On the other hand, subtle atomic 
rearrangements have been observed upon Au substitution 
for Ti2 over the planes perpendicular to [010] direction. 
In this structure, Til-01 bonds are still unequal (2.13 
and 1.82 A) making a Til-Ol-Til angle of 149.3°. Au 
substitution for a third layer Ti cation, namely s3-2x2, 
dislocates the nearby 03s in the same way as s2-2x2 
does. This suggests that Au-Ti replacement for deeper 
lying cations will cause similar local reconstruction. 

The electronic structure of s2-2x2 model in Fig.[^Ks2) 
have similar characteristics with that of sl-2x2. Both 
of them has an unoccupied defect state arising chiefly 
from Au impurity. However, the one for the s2 case lies 
higher and so partly inside CB at around T. It is slightly 
dispersed along JK and JT causing the band to be a 
quasi-lD band as opposed to the quasi-2D band of si 
case. This also shows itself in their LDOS structures. 
An 02 driven half-filled defect state having a maximum 
at r has been hived off from the lower lying group of 
surface-like bands. Therefore, Fermi level is at 0.58 eV 
with respect to bulk valence band top. This is related to 
the surface oxygen, 02, that is bonded to Au. The next 
band down is also of 02 character, however, it gets addi- 
tional contributions from the surface oxygen, 02, raised 
to the level of Ols, as well. In comparison with those of 
sl-2x2, the stack of surface bands disperse over a rela- 
tively higher energy range of 0.47 eV just above VBM. 
They are due mainly to axial Au-02 interaction. LDOS 
also shows contributions to bands slightly below VBM 
primarily from surface oxygens. As in the lxl case we 
have examined the deeper Au substitutions for up to bulk 
Ti which lead the half-filled defect state to move about 
1 eV upward closer to the empty band in the gap. 



5. Interstitials at 2x2 

Au can be encapsulated at the interstitial cavities in 
the subsurface layers starting from the fourth atomic 
layer for 2x2 surface as opposed to lxl. The first three 
possible structures are presented in Fig. [U(c)-(e). As 
common to all of them, Au relaxes into the midpoint 
between the two level oxygens where it also establishes 
equidistant coordination with four nearest neighbor 'Ti's 
that lie at the closest Ti atomic plane. Au interstitial 



causes local disturbances such that two second nearest- 
neighbor oxygens at the preceding or succeeding O lay- 
ers get slightly repelled out from their bulk lattice posi- 
tions due to the induced stress incorporated by the excess 
charge at the impurity site. Au implantation to the inter- 
nal cavities can not help lifting the symmetry breaking 
in the Til-01 bond distances that exists in the case of 
the bare anatase surface. 

The interstitial Au impurities in TiO2(001) surface, 
corresponding to Fig. HJc)-(e), show similar electronic 
characteristics since they represent equivalent local envi- 
ronments. For all of them, Fermi level falls in CB due 
to an impurity driven state that lies in the energy gap 
region along JKJ' and partly along the KT segment. The 
fully occupied quasi-lD defect state just below 1 eV dis- 
perses in accordance with the spatial alignment of Au-0 
coordination. Therefore, its maximum at J' for c-2x2 
and at J for d-2x2 and e-2x2, (doubly) alternates with 
Au-0 bond orientation. Lower lying nearly flat state in 
case-d is due primarily to 02 oxygens at the surface. For 
the cases c and especially e, the stack of bands just above 
VBM are very similar to those of the clean surface (see 
Fig. 2(a) in Rcf. [18]), because the region of disturbance 
due to the interstitial Au is away from the first Ti02 layer 
for these cases, and for cases with deeper interstitials. 



IV. BADER ANALYSIS 

We analyzed the electronic charge density using the 
atom in molecule (AIM) theory with a grid-based algo- 
rithm^. Computational Badcr charge results for near 
surface Ti and O atoms, and also for the Au impurity, 
are presented in TableHTlfor the combined systems as well 
as for the clean surface. For Au-TiO2(001) structures, 
atomic charge states are provided for O and Ti ions that 
are closest to gold site in order to better describe the 
local disturbance of the electronic density, posed by the 
impurity. 

The non-ionic nature of Ti-0 bonding in Ti02 sto- 
ichiomctry, due to below-nominal charge states of Ti 
cations and O anions, has been previously show n 29 ' 31 . 
Bulk ions accumulate Qo = — 1.33e and Qti = +2.66e 
Bader charges while these values are slightly lower for the 
corresponding surface layer species, being Q01 = — 1.26e 
and Qxii = +2.61e, due to bulk termination. 

Bader charge results presented in Table ILT1 for the impu- 
rity neighboring ions demonstrate the local disturbance 
of the electronic density. The fluctuation in the charge 
density around the cation sites remains minimal upon 
Au incorporation. For interstitial and adsorptional cases, 
gold seems to weakly bound to the lattice because of a 
limited charge transfer from Au to nearby oxygens. We 
obtained relatively smaller deviations in the charge states 
of deeper lying atoms implying a weaker contribution 
to surface electronic properties. When gold is substi- 
tuted for Ti cations it induces deviation in the electronic 
density also around the second nearest neighbor oxygens 
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owing to its spatially wider spread wavefunction. For in- 
stance, Au-Til substitution in sl-1 x 1 charges the third 
layer oxygen to Q03 = — 1.18e that is smaller than its 
reference value of Q03 = — 1.33e at the clean surface. 

In comparison with their reference values at the clean 
surface, a lower amount of charge accumulation at the 
anions around the impurity site has been obtained for all 
Au-TiO 2 (001) structures in TableQ]] This elicits a weak 
polarization in the covalency between gold and oxygen at 
this surface. The Au-0 and Ti-0 interaction strengths 
can be compared based on the charge transfers in be- 
tween. The electron depletion from Ti to O is clearly 
much larger than that from Au to O suggesting that Ti- 
O bond is stronger. 



A. Au dimer on 2x2 

The weakness of Au-0 bond polarization, in and on 
TiC>2 surface, compared to that of Ti-0 opens the possi- 
bility of an important interaction between two adjacent 
Au atoms. Clearly, Au driven impurity bands disperse 
relatively stronger when they are incorporated to 1 x 1 
surface unit cell in Fig [1] Therefore, we considered all 
possible adsorption configurations of Au dimer on 2x2 
surface in order to investigate the effect of Au-Au inter- 
action on the electronic structure of the combined sys- 
tem. Fig. [3] shows the minimum energy Au dimer struc- 
ture on the anatase TiO2(001)-2x2 surface. This is ob- 
tained from an initial configuration in which one Au atom 
is placed over the Til-Ol bond in bridge position and the 
other one is located over the next Til-Ol bond on the 
back Ti-O-Ti row. They are attracted to each other due 
to Au-Au interaction. Relatively speaking with respect 
to the structure shown in left-top of Fig [3l as they come 



TABLE II: Bader charge analysis of Au-TiC>2(001) anatase 
systems. Atom labels follow Fig. [TJ Fig. and Fig. [3] The 
values represent the valence charge states for lattice atoms 
that are closest to Au impurity site. 



Model 


Ol 


02 


03 


Til 


Ti2 


Au 


clean 


-1.26 


-1.34 


-1.33 


+2.61 


+2.64 




a-lx 1 


-1.20 


-1.34 


-1.33 


+2.58 


+2.64 


-0.03 


sl-lxl 


-0.70 


-0.85 


-1.18 




+2.64 


+1.38 


s2-l x 1 


-1.24 


-1.17 


-0.94 


+2.58 




+1.54 


a-2x2 


-1.15 


-1.34 


-1.32 


+2.55 


+2.65 


-0.04 


b-2x2 


-1.18 


-1.35 


-1.33 


+2.57 


+2.62 


+0.34 


c-2x2 


-1.26 


-1.33 


-1.24 


+2.55 


+2.56 


+0.34 


d-2x2 


-1.25 


-1.34 


-1.30 


+2.60 


+2.56 


+0.34 


e-2x2 


-1.26 


-1.34 


-1.34 


+2.60 


+2.58 


+0.38 


sl-2x2 


-0.95 


-1.11 


-1.26 




+2.63 


+1.39 


s2-2x2 


-1.20 


-1.03 


-1.27 


+2.63 




+1.32 


aa-2x2 (Aul) 


-1.23 


-1.35 


-1.32 


+2.56 


+2.64 


-0.07 


aa-2x2 (Au2) 


-1.12 


-1.35 


-1.33 


+2.54 


+2.64 


+0.05 




[010] r j k r r k 



FIG. 3: (color online) Au dimer on anatase TiO2(001)-2x2 
surface (aa-2x2). Side views of the impurity embedded slab 
model are shown in the left panel where all distances are given 
111 A. Top view of gold dimers on this surface as well as the 
corresponding electronic structure are presented on the right 
panel. 

closer reducing the dimer length to 2.55 A the one at 
the back row (Au2) pulls the nearest neighbor Ol off its 
lattice position up by 0.91 A. However, the Ol in inter- 
action with the Aul at the first row is elevated by only 
0.06 A. Surprisingly, any geometry that the surface oxy- 
gens bonded symmetrically to the dimer is energetically 
unfavorable. 

We also examined the possibility of a dissociative ad- 
sorption for various configurations which corresponds to 
0.5 ML coverage on this surface. Our results suggest that 
such a dimer dissociation on the anatase surface is ener- 
getically not preferable signifying the strength of Au-Au 
attraction compared to Au-0 and Au-Ti interactions. 

Electronically, gold dimer adsorption at 2x2 unit cell 
drives the anatase (001) surface into a semiconducting 
state such that it reduces the band gap by 1.04 eV rela- 
tive to anatase bulk VBM. The impurity band just below 
the Fermi energy has a minimum at T and a maximum 
at J' which gives rise to an indirect gap. This otherwise 
almost fiat going state lies higher in energy from a group 
of five surface-like bands with a separation of 0.34 eV 
at r. It is mainly characterized by Au-Au and Aul-Ol 
interactions with a partial mixing from the unsaturated 
surface oxygens. The second defect states disperses con- 
jugately to the first one having a maximum at T and a 
minimum at J'. The surface oxygen that is bonded to 
Au2 mainly contributes to the second and the third de- 
fect states. The next impurity band that has a minimum 
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at J and a maximum at J' originates from Au2-01 bond- 
ing. The LDOS of unsaturated 01 ions extend over the 
VB, also giving partial contribution to all of the defect 
states. 



V. THERMODYNAMIC STABILITY 

We analyzed the thermodynamic stability of gold in- 
corporated T1O2 systems, considered in this work, as- 
suming Au and Ti exchange between the substrate and 
the surrounding gas phase to account for the formation 
of substitutional cases. 

We assume experimental conditions so that Au is cho- 
sen to be in equilibrium with its metallic bulk phase and 
redefine the zero point of energy at the maximal value of 
/JTi by introducing A^Ti = MTi — A*Ti ■ The relative sur- 
face GFEs, A7, (or the formation energies) are plotted 
as a function of A^xi in Fig. [4] over its allowed range of 
values for 12 surface structures considered in this study. 
The lower limit for Ti chemical potential has been ob- 
tained to be —8.90 eV due to the fact that the surface 
must keep in thermal equilibrium with the bulk anatase 
TiC>2 phase. 

Our prescription of stability of the phases can pre- 
dict the structural reconstructions rather than thermally 
stimulated formations due to the approximation made by 
omitting the entropy terms in the surface GFE. 

The results presented in Fig. [4] suggest that Au incor- 
poration at lxl surface is not favorable. For instance, 
Au substitution for the second layer Ti at lxl anatase 
(001) surface (s2-lxl) is energetically the least prefer- 
able phase among the structures considered. Similarly, 
sl-lxl case is also thermodynamically unstable. Their 
formation energies further increase linearly toward Ti— 
rich experimental conditions. Au adsorption for 1 ML 
coverage would require an energy of 1.23 eV/lxl rela- 
tive to that of the clean surface. 

Au intcrstitials in the 2x2 surface appears to be un- 
favorable compared to bulk terminated bare anatase 
TiO2(001) by the formation energies of 0.44, 0.35 and 
0.35 eV/lxl for c, d, and e cases, respectively Surpris- 
ingly, Au implantation into the interstitial cavities deeper 
than the fourth atomic layer is slightly more preferable 
at this surface. These calculated formation energies for 
Au interstitials in anatase surface are in good agree- 
ment with the experimental results of Perkas et alJ^. 
In an attempt to explain Au induced crystallization (IC) 
of anatase through gold insertion into Ti02 which re- 
quires thermal treatment at about 80°C by sonication, 
they stated that Au diffuses into the support layers and 
forms and intermixed layer through a multiphase diffu- 
sion across the gold support interface. 

Substitutional gold impurities tend to be linearly more 
unstable toward Ti-rich conditions. It is also reasonable 
to see that Au substitution for a second Ti layer cation 
is energetically more unstable than that for the surface 
layer Ti over the full range of allowed /ixi values. On the 
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■/ 
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-8.9 -6.0 -3.0 0.0 
O-rich -+ A^iTi(cV) *- Ti-rich 

FIG. 4: (color online) Normalized relative surface Gibbs free 
energies of Au-TiO 2 (001)-lx 1 and Au-TiO 2 (001)-2x2 struc- 
tures as a function of and over the allowed range of Ti chem- 
ical potential. Au impurity ions are chosen to be in thermo- 
dynamic equilibrium with ccp Au bulk phase. 



other hand, Au(Ti) substitution has been expected to be 
more preferable in a Ti-poor environment. Indeed, both 
of the substitutional cases become more stable than the 
clean relaxed surface under O-rich conditions when the 
substrate establishes thermodynamic equilibrium with 
the surrounding oxygen gas phase. Fig. [4] shows that 
such a surface structure can be realized under the limits, 
A/zxi < -7.84 cV and A^ T i < -7.14 cV for s2 and si 
cases, respectively. Furthermore, Au(Ti) substitution for 
the surface Ti atom proves to be the most stable phase 
for —8.90 < A/iTi < —7.98 eV that corresponds to a 
Ti-poor environment. 

1ML gold adsorption is not preferable with a forma- 
tion energy of 1.23 eV relative to that of the clean sur- 
face. On the other hand, Au adsorbates at 2x2 unit 
cell as well as the dimer structure in Fig. [3] happen to 
be even more stable than the formation of bulk termi- 
nated bare (001) surface which is known to be under 
tensile stress^ whereby Au bonding causes a lowering 
of the surface energy by releasing the surface strain^ 
We calculated the corresponding relative surface GFEs 
as -0.17, -0.21 and -0.18 eV/lxl for a-2x2, b-2x2 
and aa-2x2 structures, respectively. Therefore, highly 
symmetrical structure in Fig. [2](b) appears to be the most 
stable phase over a wide range of experimental situations 
within —7.84 < A/ixi < from Ti-low to Ti-rich equi- 
librium conditions. However, the differences between the 
computed GFE values of these three adsorption modes 
are, surprisingly, not significantly large and can be ac- 
cepted within the limits of a computational accuracy. In 
this sense, one cannot conclude that a-2x2 and aa-2x2 



9 



are considerably less stable than b-2x2 is. The results 
rather suggest that these three phases might co-exist in 
an experimental situation where, interestingly, single Au 
atom adsorption is metallic while dimcr structure is semi- 
conducting. Therefore, overall behavior is expected to be 
conducting in nature. 

VI. CONCLUSIONS 

In conclusion, we studied the role of Au mediated lat- 
tice relaxations as well as the effect of the impurity itself 
on the electronic structure of anatase TiO2(001) for lxl 
and 2x2 surface unit cells. We, additionally, considered 
gold dimcr adsorption on TiO2(001)-2x2 to account for 
the strength of Au-Au interaction over Au-0 attraction. 

In comparison with Ti, due to its spatially wider spread 
wavefunction, gold maintains a relatively larger distance 
to the oxygens with which it interacts. In addition, Au 
exhibits an ample coordination with the neighboring oxy- 
gens and, therefore, disturbs local lattice structure. Elec- 
tronically, an Au adsorbate transfers a limited amount of 
charge to the nearest neighbor oxygens that reduces the 
ionicity leading to a relatively weaker impurity-oxygen 
interaction. On the other hand, Au can attain higher 
valence charge states in the case of interstitial and sub- 
stitutional incorporations due to increased coordination 
number with the neighboring oxygens. For 1ML impurity 
concentration, gold interstitial enforced local distortion 
causes significant reconstruction in the form of interlayer 
segregation due to shortened Au-Au interaction distance. 
These impurity mediated disturbances derive new defect 
states and strongly disperse the existing surface bands. 

Single Au impurity brings a half-filled impurity state 
into the band gap of TiC>2(001) which pins the Fermi 
level leading to metallization, in addition to other filled 
surface and impurity driven bands within the gap. In 
the adsorptional and interstitial cases, this state is de- 



rived from the CB. For Au interstitials as well as for 
b-2x2, Fermi energy falls inside the CB. However, sub- 
stitutional Au makes it up close to the VB also deriving 
empty impurity states higher in energies toward the CB. 
The dispersion of the defect states depend on the im- 
purity concentration and the Au-0 interaction strength. 
Therefore, lxl cases exhibit the strongest dispersions, 
and the adsorptional Au incorporations at 2x2 surface 
lead to flat-like defect states which represent the weakest 
Au-0 interactions. 

All Au-incorporated TiO2(001) structures at lxl sur- 
face unit cells are unstable. Under strongly oxidizing, Ti— 
poor conditions gold tends to substitute surface Ti ions 
at 1/4 ML coverage. On the other hand, for a wide range 
of experimental situations, from moderate to strongly re- 
ducing conditions, Au adsorption is thermodynamically 
more stable at the 2x2 support surface. Single Au adsor- 
bates at this coverage bring out new acceptor sites as they 
reduce the surface leading to metallization. Therefore, 
an additional gold adsorption prefers the dimer structure 
rather than being dissociative. Gold dimcr supported by 
the slab at 2x2 surface is almost as stable as the sin- 
gle Au adsorption signifying the strength and the role of 
Au-Au interaction which drives the system into semicon- 
ducting regime. Consequently, in an experimental situa- 
tion which realizes the co-existence of these two phases, 
the overall behavior is expected to be conducting. 
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